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Abstract 



A search for microscopic black holes in pp collisions at a center-of-mass energy of 
7TeV is presented. The data sample corresponds to an integrated luminosity of 
4.7fb^^ recorded by the CMS experiment at the LHC in 2011. Events with large 
total transverse energy have been analyzed for the presence of multiple energetic 
jets, leptons, and photons, which are typical signals of evaporating semiclassical and 
quantum black holes, and string balls. Agreement with the expected standard model 
backgrounds, which are dominated by QCD multijet production, has been observed 
for various combined multiplicities of jets and other reconstructed objects in the final 
state. Model-independent limits are set on new physics processes producing high- 
multiplicity, energetic final states. In addition, new model-specific indicative limits 
are set excluding semiclassical and quantum black holes with masses below 3.8 to 
5.3 TeV and string balls with masses below 4.6 to 4.8 TeV. The analysis has a substan- 
tially increased sensitivity compared to previous searches. 
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1 Introduction 

One of the most spectacular predictions of theories with low-scale quantum gravity is the pos- 
sibility of microscopic black hole (BH) production in proton-proton collisions at the high en- 
ergies offered by the Large Hadron Collider (LHC) [1, 2j. Such models are motivated mainly 
by the puzzling large difference between the electroweak scale (~0.1 TeV) and the Planck scale 
(Mpi ~ 10^^ TeV), known as the hierarchy problem. In this analysis, we focus on black hole 
production in a model with n large, flat, extra spatial dimensions (ADD model) HsHH. In this 
and in other models, the fundamental scale of new physics in n extra dimensions is given in 
terms of a multidimensional Planck scale Md, such that Mj^^^ oc Mp^R^", where R is the size of 
extra dimensions. Some of the conclusions also apply to black holes in the Randall-Sundrum 
model ISJEI/ with a single warped extra dimension. 

This analysis extends a previous search |j7| for short-lived microscopic black holes carried out 
by the Compact Muon Solenoid (CMS) Collaboration in 2010. The present search is based on 
the full 2011 data sample, which corresponds to an integrated luminosity of 4.7 ± 0.2 fb^^ [8, 9J 
at a center-of-mass energy of 7 TeV. The details of the analysis method and the underlying 
theory, as well as the detailed description of black hole evaporation models, can be found in the 
original publication [7J. Typically, microscopic black holes are characterized by the production 
of a large number of energetic final-state particles, 75% of which are jets. Searches for black 
holes have also been performed by the ATLAS Collaboration IITOl[TT| . 

We present our results in terms of model-independent limits on the cross section times the 
branching fraction into a multiparticle final state. Further, we interpret the results in terms of 
a set of benchmark black hole models. The analysis also extends the previous CMS search for 
semiclassical black holes to probe for other quantum gravity objects such as string balls [12] 
and quantum black holes |T3l. It is commonly accepted that the semiclassical approximation 
is valid when the black hole mass is some 3-5 times larger than the Md- The string balls are 
hypothetical precursors of semiclassical black holes in an extreme quantum limit, when the 
mass of the object is close to the Planck scale. In cases where the semiclassical approximation 
no longer holds, string balls may offer a more realistic description of black hole formation 
and decay. String balls are described by the string scale Ms and string coupling gs- These 
objects would evaporate similarly to black holes, except that the evaporation would occur at 
the Hagedorn temperature, which does not depend on the string-ball mass ||T2l [14], unlike 
the Hawking temperature [15J, which decreases as the black hole mass increases. Another 
possibility is that a light black hole with mass close to the Planck scale may evaporate faster 
than it thermalizes, resulting predominantly in a nonthermal decay into a pair of jets, rather 
than into high-multiplicity final states [ISKISKlZ]- We search for production of these objects, 
referred to as quantum black holes, and for their decay in both the ADD scenario and in the 
Randall-Sundrum model of low-scale gravity with a single (n = 1) compact extra dimension. 

2 The CMS detector 

A detailed description of the CMS experiment can be found elsewhere |[T8|| . The CMS detector 
consists of a 3.8 T superconducting solenoid enclosing a silicon pixel and strip tracker, a crystal 
electromagnetic calorimeter (ECAL), and a brass-scintillator hadronic calorimeter (HCAL). The 
finely segmented ECAL employs lead-tungstate crystals with transverse dimensions: A;/ x 
= 0.0174 X 0.0174. The HCAL cells are grouped in projective towers, of granularity A?/ x 
A(^ = 0.087 X 0.087 at central rapidities and increasing progressively in the forward region. 
Here, (p and 9 are the azimuthal and polar angles, with 6 measured with respect to the direction 



2 



3 Event reconstruction and Monte Carlo samples 



of the counterclockwise beam. The pseudorapidity f] is defined as — ln[tan(0/2)]. Muons are 
measured in the pseudorapidity range \rj\ < 2.4 in gas-ionization detectors embedded in the 
steel return yoke. 

The CMS trigger system, used to select the most interesting events, consists of two levels. A 
first level (LI), composed of custom electronics, uses information from the calorimeters and 
muon detectors to decrease the event rate to 80 kHz. A software-based High Level Trigger 
(HLT) further decreases the event rate to 350-400 Hz for data storage. 

As in the previous analysis, we use data collected with a suite of Hj triggers, where Hj is 
defined as a scalar sum of the transverse energies (Ej) of the jets above a thresholdj^ There 
have been changes introduced in the trigger logic both at LI and HLT since 2010. We now use 
jets corrected for the calorimeter response to calculate the Hj variable at the HLT (uncorrected 
jets are still used at LI). Also, the minimum Hj thresholds, as well as the minimum jet Ej 
requirement for a jet to be counted towards Hj, have been increased to account for pileup 
effects and to allow for the increased instantaneous luminosity of the LHC. This minimum 
jet Et threshold is lOGeV at LI, and 40GeV at the HLT. The minimum Hj threshold at the 
HLT varies between 350 and 650 GeV, depending on the instantaneous luminosity. Only jets 
reconstructed at central pseudorapidities \rj\ < 3.0 are used in the Hj calculations at LI and at 
the HLT for the 2011 data-taking period. The trigger is fully efficient for events with Hj above 
800 GeV. In order to explore all possible black hole decay modes, the entire analysis was also 
repeated using data collected with multimuon or missing transverse energy (£^^^*) triggers, 
but this yielded no events consistent with the expected black hole production. 

3 Event reconstruction and Monte Carlo samples 

Jets are reconstructed offline using energy deposits in the HCAL and ECAL, which are clus- 
tered using an infrared-safe anti-A;^ algorithm with a distance parameter of 0.5 FT9l-l2T|. Qual- 
ity criteria are applied to jets in order to remove calorimeter noise and noncollision back- 
ground ||22|. We require jets to have Ej > 20 GeV and to be reconstructed within \rj\ < 2.6. 
Further, jet energies are corrected for calorimeter nonuniform response with correction factors 
derived from Monte Carlo (MC) simulation and dijet events from collision data [22J. The trans- 
verse energy resolution for jets AEj/ Ej is better than 15% in the range considered. We recon- 
struct E^'^^ as the modulus of the negative vector sum of transverse energies in the individual 
calorimeter towers and is further corrected for the jet energy scale and for muon momenta 
measured in the trackers f23l. 

Electrons and photons are reconstructed as isolated energy deposits in the ECAL with a shape 
consistent with that expected for electromagnetic showers. Electrons are required to have a 
track matched to the calorimeter cluster, while photons are required to have no matching hits 
in the silicon layers. Electrons and photons are selected with Ej > 20 GeV and are required to 
be reconstructed in the fiducial barrel {\t]\ < 1.44) or endcap (1.56 < \f]\ < 2.4) regions. The 
ECAL has excellent energy resolution, for example contributing 1 GeV to the observed width of 
the Z boson using e^e^ pairs with low bremsstrahlung loss, measured in the barrel region [24j. 

Muons are reconstructed as matched tracks in the muon spectrometer and the silicon tracker. 
Muons are selected with \rj\ < 2.1 and pj > 20 GeV, and are required to be isolated from 
other tracks. Requiring the muons to have distance of the closest approach < 0.2 cm helps to 
suppress backgrounds from cosmic-ray muons. Here, the distance of the closest approach is 
defined as the shortest distance between the beam line and the direction of an object in the 



Energetic electrons and photons are also reconstructed as jets at the trigger level. 
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transverse plane. Performing a combined fit to the track segments measured in the silicon 
tracker and the muon system results in a transverse momentum resolution between 1% and 
5% for pj values up to 1 TeV. The minimum separation between any two objects (jet, lepton, or 
photon) is required to be AR = ^/{Acp)^ + (Af/)^ > 0.5. 

Simulated samples of semiclassical black hole events are produced using a parton-level Black- 
Max II25I (v2.01) and CHARYBDIS (IIHSZI (CHARYBDIS 2, vl.0.3) MC generators, followed by 
the parton showering simulation with PYTHIA |28| (v6.420), and a fast parametric simulation 
of the CMS detector |29|. The MSTW20081o68cl [30| parton distribution functions (PDF) are 
used for generating the signal samples. The BlackMax and CHARYBDIS generators calculate 
the total cross section a from geometric considerations, assuming that a cx yrrg, where rg is 
the Schwarzschild radius [1. 2J. The BlackMax generator uses the approximation of rotating 
black holes, which includes an additional factor that depends on the number of extra dimen- 
sions n. The CHARYBDIS generator incorporates a more detailed model based on Yoshino- 
Rychkov corrections |l3Tl|32l to the pure geometrical cross section, resulting in production cross 
sections that are lower by a factor of 1.36, 1.59, and 1.78 compared to those from BlackMax 
for n = 2, 4, and 6, respectively. These scale factors can be used to interpret the results obtained 
with one framework in terms of the other. 

Certain models are supported by both generators: rotating and nonrotating black holes, and 
black holes with mass and angular momentum loss prior to evaporation. This loss is set to 
be 10% in the BlackMax generator, while it is estimated in the Yoshino-Rychkov model and 
varies from 18% to 30% for n = 2 to 6 in the CHARYBDIS generator. In addition, we use CHARYB- 
DIS to simulate black hole evaporation resulting in a stable remnant with mass equal to the 
multidimensional Planck scale Md, or a boiling remnant (unique to the CHARYBDIS generator). 
Both of these scenarios represent alternative descriptions of the final stage of the black hole 
evolution. In the case of a stable remnant, the terminal stage of a black hole is a noninteracting 
remnant with a mass of order Md', in the case of a boiling remnant, a black hole undergoes a 
transformation into a string ball at a mass close to Md with subsequent evaporation at a fixed 
temperature. We produce a number of string-ball samples using the BlackMax generator. Fi- 
nally, the QBH (version 1.03) matrix-element generator f33| with CTEQ6L PDF set [34 1 is used, 
followed by the parton showering simulation with PYTHIA and fast simulation of the CMS de- 
tector, to produce quantum black hole samples. Table [l] summarizes the models used in this 
search. 



Table 1: Signal Monte Carlo samples and generators used in the analysis. 



Sample description 


BlackMax 


CHARYBDIS 


QBH 


nonrotating BH 


YES 


YES 


NO 


Rotating BH 


YES 


YES 


NO 


Rotating BH with mass and 








angular momentum loss 


YES (10% loss) 


YES (18 - 30% loss) 


NO 


Rotating BH, low multiplicity regime 


NO 


YES 


NO 


Boiling remnant 


NO 


YES 


NO 


Stable remnant 


NO 


YES 


NO 


String balls 


YES 


NO 


NO 


Quantum BH 


NO 


NO 


YES 
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4 Analysis method 



4 Analysis method 

The total transverse er\ergy is used to separate black hole candidate events from backgrounds. 
A variable Sj is defined as a scalar sum of the transverse energy (Ej) of individual objects: 
jets, electrons, photons, and muons passing the selections described above. Only objects with 
Ej > 50 GeV enter into the sum for the calculation of Sj and count towards the final-state 
multiplicity N. This rather high minimum transverse energy requirement makes the analysis 
insensitive to the jets from pileup and reduces the SM backgrounds by a few orders of magni- 
tude, while being fully efficient for black hole decays. We further add the measured E™^^ in 
the event to the Sj, if E^^^ > 50 GeV. Generalization of the Sj definition to include £^'®® is 
important for testing black hole models with a significant amount of missing energy such as 
models with a stable noninteracting remnant. A spurious E^^^ may arise in an event as the 
result of mismeasurement of the jets. However, we have checked that the consequent effect on 
St of double counting energy in both the jet and £^'®® contributions is negligible. Note that by 
construction, particle misidentification does not affect the total transverse energy in the event 
considerably. 

Depending on the details of the black hole evaporation, a large variation of particle multiplicity 
in the final state, and a large range of missing transverse energies are possible. While resulting 
in quite different signatures, these variations typically have very little effect on the value of 
St in the event. A recent work on quantum-gravity black holes | |35| also suggests a larger 
number of softer particles produced in black hole evaporation than in the semiclassical case, 
further emphasizing the importance of St as a largely model-independent variable for black 
hole searches. 

The main background to black hole production arises from QCD multijet production, which 
dominates the event rates at large St. Smaller backgrounds come from 7/W/Z+jets and tt 
production. These smaller backgrounds are negligible at large values of St and contribute less 
than 1% to the total background after the final selection. We estimate their contribution from 
MC simulation, using the MadGraph [36J leading-order parton-level event generator (with 
up to three extra partons included in the simulation) with the CTEQ6L PDF set followed by 
PYTHIA 1281 parton showering and full CMS detector simulation via GEANT4 [37J. For the 
dominant QCD background, however, we estimate backgrounds from the observed data using 
the St multiplicity invariance method [7J. This method relies on the independence of the shape 
of the St spectrum on the number of final-state objects N; an empirical observation extensively 
checked by using various MC samples (alpgen and pythia) as well as low-multiplicity data. 
The origin of this invariance lies in the collinear nature of the final-state radiation, which typi- 
cally does not change the total transverse energy in the event; hence the independence of the St 
spectrum of the jet multiplicity for the QCD background. This invariance allows us to predict 
the shape of the St spectrum for any number of objects using the dijet data, which has been 
studied extensively for presence of new physics in dedicated analyses I|38ti40| . 

We use low-multiplicity data with N = 2 and N = 3 to obtain the background shape by fitting 
the St distributions between 1200 and 2800 GeV with the ansatz function Po{l + x)^^ /x^^+^s 
which is shown in Fig. [l] as a solid line. No evidence of new physics has been observed in 
this region in a dedicated analysis Il38l . To estimate the systematic uncertainty of the method, 
the same St distributions are fitted with two additional functions, Pq/ (^"1 + Pi^ + x^)^^ and 
Pq/ {Pi + x)^^. Thus, an envelope of functions is formed (shown as the shaded area in Fig. [ij 
and is used as the systematic uncertainty. 

Figure |2] shows the fit result of the background prediction for the inclusive samples with high 
object multiplicity events. Here, the shape of the St distribution obtained from the N = 2 
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Figure 1: Distribution of the total transverse energy, Sj, for low-multiplicity events with mul- 
tiplicity: a) N = 2 and b) N = 3 photons, electrons, muons, or jets in the final state. Ob- 
served data are depicted as points with error bars; solid line with a shaded band is the back- 
ground prediction and its systematic uncertainty. Non-QCD backgrounds are shown as filled 
histograms (not stacked). Also shown is the black hole signal for three parameter sets of the 
BlackMax nonrotating black hole model, demonstrating that signal contamination in the fit 
region of 1200 - 2800 GeV would be small. 

sample is normalized to the observed data in the range 1800 to 2200 GeV, where no signal 
contribution is expected. Also shown are the expected semiclassical black hole signals for three 
parameter sets of the BlackMax nonrotating black hole model. The results are presented 
separately for six different values of the minimum final state multiplicity. The data agree with 
the background shapes from the low-multiplicity samples and do not exhibit evidence for new 
physics. Figure |3] shows a similar comparison of the experimental Sj distribution with the 
predicted signal for three parameter sets of the QBH quantum black hole model. In this case 
the comparison is shown separately for just two values of the minimum final state multiplicity, 
reflecting the different decay characteristics expected for quantum black holes compared to 
semiclassical black holes. 



5 Results 

In order to set exclusion limits on black hole production, we assign systematic uncertainties on 
the background estimate varying from 3% to 300% in the Sj range used in this search. These 
uncertainties are dominated by the uncertainties from using various fit ansatz functions (2%- 
300%), which are added in quadrature to the second-largest contribution, which arise from the 
normalization statistical uncertainty (2%-21%). The integrated luminosity is measured with 
4.5% uncertainty [8. 9J utilizing information from the forward calorimeters. The signal uncer- 
tainty is dominated by the jet energy scale uncertainty of ~ 2% [22J, which translates into 2% 
uncertainty on the signal. An additional 2% uncertainty on the signal acceptance comes from 
the variation of acceptance obtained with the default MSTW20081o68cl PDF library and PDFs 
within the CTEQ61 and CTEQ66 error sets L34J . 

Given the significant model dependence of the black hole production cross section and decay 
patterns, it is not practical to test all different variations of model parameters, offered by recent 
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Figure 2: Distribution of the total transverse energy, Sj, for events with multiplicity: a) N > 3, 
b) N > 4, c) N > 5, d) N > 6, e) N > 7, and f) N > 8 objects (photons, electrons, muons, or 
jets) in the final state. Observed data are depicted as points with error bars; the solid line with 
a shaded band is the background prediction and its systematic uncertainty. Also shown are the 
expected semiclassical black hole signals for three parameter sets of the BlackMax nonrotat- 
ing black hole model. Here, is the minimum black hole mass, Md is the multidimensional 
Planck scale, and n is the number of extra dimensions. 
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Figure 3: Distribution of the total trar\sverse energy, Sj, for events with multiplicity: a) N > 2 
and b) N > 3 objects in the final state. Observed data are depicted as points with error bars; 
the solid line with a shaded band is the background prediction and its systematic uncertainty. 
Also shown are the expected quantum black hole signals for three parameter sets. Here, 
is the minimum quantum black hole mass, Md is the multidimensional Planck scale, and n is 
the number of extra dimensions. 



black hole event generators, in a dedicated search. This study considers some 700 different 
signal MC samples, yet it does not come close to spanning the entire parameter space of the 
models; scaling the number of signal samples up and presenting the results for every model 
therefore becomes impractical. Hence, we first present the results of our search in a generic, 
model-independent way, which would allow others to probe additional models using parton- 
level MC information, possibly augmented with a very basic detector simulation. To facilitate 
such an approach, we provide model-independent limits on the cross section times the accep- 
tance for new physics production in high-Sj inclusive final states for N > 3,4,5,6,7, and 8. 
The limits are set using a modified frequentist CLg method ||4TI|42| with log-normal prior used 
to marginalize nuisance parameters in the likelihood function. 

Figure |4] shows 95% confidence level (CL) limits from a counting experiment placed on the 
experimentally reconstructed value Sj > S^"^ as a function of S^'", which can be used to test 
models of new physics that result in these final states, including (but not limited to) an even 
broader variety of black hole models than we covered in this analysis. The 95% CL limits from 
2010 data [7] are also shown in Fig.|4]for comparison. The present model-independent limits are 
roughly 0.6 fb for high values of Sj, representing a two orders of magnitude improvement over 
the limits reported in our first publication [7J. Given the higher statistics of the 2011 sample, 
we are also able to extend these limits to the N > 6,7, and 8 cases. 

For a specific subset of the black hole models fZl that are being probed, we also set dedicated 
limits on semiclassical and quantum black hole and string-ball production performing count- 
ing experiments using optimized Sj and N selections. It should be noted that the semiclassical 
approximation used for deriving the cross section within respective benchmark scenarios is 
expected to break down for many of the points probed, a point emphasized in a recent cri- 
tique [43J. Thus, these limits should be treated as indicative, rather than precise. 

The signal (S) significance is optimized in the presence of background (B) using a test statis- 
tic S/ \/S + B for each set of model parameters. The optimum choices of Sj and N for a few 
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Figure 4: Model-independent 95% CL upper cross section limits for counting experiments with 
St > S^"" as a function of S^^^ for events with multiplicity: a) N > 3, b) N > 4, c) N > 5, d) 
N > 6, e) N > 7, and f) N > 8. The blue solid (red dotted) lines correspond to an observed 
(expected) limit for nominal signal acceptance uncertainty of 5%, compared to observed (ex- 
pected) limits obtained with 2010 CMS data and shown as blue dashed (red dash-dotted) line. 
The green and yellow bands represent one and two standard deviations from the expected 
limits. 



9 



illustrative benchmark scenarios are listed in Table |2| as well as the predicted number of back- 
ground events, expected number of signal events, and the observed number of events in data. 
The corresponding cross section limits at the 95% CL from the counting experiments are shown 
in Fig.|5]for a small subset of characteristic signal points. 

Table 2: Details of some of the BlackMax nonrotating black hole model parameters probed 
in the analysis. Shown are: multidimensional Planck scale (Mo), minimum black hole mass 
(Mbh)/ number of extra dimensions (n), corresponding leading order cross sections {a), and 
optimal selections on the minimum decay multiplicity (N > N"^^) and minimum Sj, as well 
as signal acceptance (A), expected number of signal events (N^^S)^ number of observed events 
in data (W^^'^), expected background with its systematic uncertainty, and observed ((7^^) 

and expected ((cr^^)) limit on the signal at 95% CL. 



Md 


Mbh 


n 


a 






A 








(7^5 




(TeV) 


(TeV) 




(Pb) 




(TeV) 


(%) 








(pb) 


(pb) 


1.5 


3.0 


6 


26 


3 


1.9 


88.0 


110000 


5999 


5970 ± 180 


0.092 


0.091 


1.5 


3.5 


6 


5.0 


3 


2.2 


88.0 


21000 
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By translating the cross section limits into ADD model expectations, we can exclude the pro- 
duction of semiclassical black holes with a minimum mass varying from 3.9 to 5.3 TeV for values 
of the multidimensional Planck scale Md < 4 TeV and a number of extra dimensions n < 6 at 
95% CL (see Fig. |6]|. The excluded minimum masses of quantum black holes are in the 3.8- 
5.2 TeV range for Md up to 4 TeV and are shown in Fig. [7| The 95% CL limits on the string-ball 
production cross section as a function of the minimum mass of the string ball are shown in 
Fig. |8] We exclude string-ball production with a minimum mass from 4.6 to 4.8 TeV, depend- 
ing on the model. Despite the caveats mentioned above, we consider it useful to present these 
results, since at present there are no alternative quantitative calculations in the regime where 
the semiclassical approximation breaks down. Furthermore, although the predicted cross sec- 
tion is very sensitive to changes in the model of black hole production and decay, varying the 
model assumptions results in only moderate changes to the mass limit, because of the exponen- 
tial dependence of the cross section on the black hole mass. Nontheless, we emphasize that the 
model-independent limits set in this Letter should be used in the regime when the semiclassical 
approximation fails in order to obtain more reliable predictions. 

6 Conclusions 

An update has been presented of an earlier dedicated search for black holes at the LHC JZl 
and new model-independent limits have been set on the production of energetic multiparticle 
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Figure 5: Cross section limits at 95% CL from the counting experiments optimized for various 
black hole parameter sets (solid lines) compared with signal production cross sections from the 
BlackM AX generator (dashed lines) as a function of minimum black hole mass. 

final states, which can be used to constrain a large variety of models of new physics. For the 
benchmark models, the following can be concluded from Fig.|6j (i) since the excluded masses of 
rotating and nonrotating black holes are similar, the effect of black hole spin on the sensitivity 
of the search is small; (ii) in case of energy/ momentum loss due to gravitational radiation not 
trapped by the forming event horizon, the excluded black hole masses are ~10% lower than in 
the case of no losses; (iii) the choice of Sj as a discriminating variable makes the results largely 
insensitive to the details of the last stage of black holes evaporation, whether a stable remnant is 
formed or not. Numerically, the limits on the minimum semiclassical and quantum black hole 
and string-ball masses are in the range 3.8 to 5.3 TeV for a wide range of model parameters. 
These are the most restrictive limits on black hole production set at hadron colliders to date. 
Further extension of this search will be possible when the LHC energy is increased in the future. 
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Figure 6: Minimum black hole mass excluded at 95% CL as function of the reduced Planck 
scale for various BlackMax black hole models without the stable remnant and number of 
extra dimensions of two, four, and six (Top left). The minimum black hole mass, excluded at 
95% CL, as function of the reduced Planck scale for various CHARYBDIS black hole models with 
or without the stable remnant and number of extra dimensions n of (top right) two, (bottom 
left) four, and (bottom right) six. The areas below each curve are excluded by this search. 
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Figure 8: Cross section limits at 95% CL from the counting experiments optimized for various 
string-ball parameter sets (solid lines) compared with signal production cross section (dashed 
lines) as a function of minimum string-ball mass. Here, Md is the multidimensional Planck 
scale, Ms is the string scale, and gs is the string coupling. 
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